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Most  Corps  constructed  tunnels  are  either  water  conveyance  or  trans- 
portation tunnels  driven  through  rock.  This  paper  is  limited  to  ways  in 
which  the  computer  can  be  used  to  aid  the  structural  engineer  in  the  de- 
sign of  such  tunnels  A The  engineer  must  design  the  temporary  (steel  sets, 
rock  bolts,  spiling, riftiotcrete)  and  permanent  tunnel  supports.  The  paper 
discusses  programs  available  for  designing  steel  sets,  concrete  linings, 
and  steel  tunnel  linersu  surveys  the  finite  element  method,  and  suggests 
developing  programs  for  rock  bolts,  spiling,  and  shotcrete.  The  paper 
concludes  that  although  the  computer  cannot  be  used  to  automate  the  de- 
sign of  tunnel  supports,  it  can  be  used  to  aid  in  making  sound  engineer- 
ing judgements. 
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In  December  197*+,  the  Automatic  Data  Processing  (ADP)  Center, 
Waterways  Experiment  Station  (WES)  submitted  a proposal  to  conduct  a 
Corps-wide  Conference  on  Computer-Aided  Design  in  Structural  Engineer- 
ing to  the  Office  of  the  Chief  of  Engineers  (OCE).  OCE  approved  the 
proposal  and  efforts  were  started  in  February  1975  to  conduct  this 
Conference.  The  Conference  was  conducted  in  New  Orleans,  Louisiana, 
22-26  September  1975  and  was  attended  by  175  engineers  from  U8  Corps 
field  offices,  OCE,  Construction  Engineering  Research  Laboratory,  and 
WES. 

This  volume  contains  the  papers  from  Specialty  Session  C,  State- 
of-the-Corps  Art  on  Single-  and  Multicell  Conduits  and  Tunnels. 

Mr.  George  W.  Henson,  Structural  Engineer,  SWTED-DT,  Tulsa  Distruct, 
was  session  chairman  and  presented  a paper.  Other  papers  were  pre- 
sented by  Mr.  Carney  M.  Terzian,  Chief,  Structural  Section,  NEDED-T, 
New  England  Division,  and  Mr.  Robert  J.  Smith,  Structural  Engineer, 
DAEN-CWE-D,  OCE. 

The  Conference  was  successful  due  to  the  efforts  of  a multitude 
of  people.  The  roles  they  played  were  different  but  they  were  all 
directed  toward  making  a concept  on  "instant  dissemination"  work.  The 
Organizing  Committee  for  the  Conference  consisted  of: 

COL  G.  H.  Hilt,  WES 
Mr.  F.  R.  Brown,  WES 
Mr.  D.  L.  Neumann,  WES 
Mr.  J.  B.  Cheek,  Jr.,  WES 

Dr.  N.  Radhakrishnan , WES — Conference  Coordinator 

Mr.  W.  A.  Price,  WES 

Mr.  G.  S.  Hyde,  WES 

Mr.  D.  R.  Dressier,  LMVD 

Mr.  W.  B.  Dodd,  LMNDE 

Ms.  E.  Smith,  LMNDE 

Mr.  L.  H.  Manson,  LMNDE 
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An  OCE  Coordinating  Committee  also  worked  enthusiastically  to 
ensure  the  success  of  the  Conference.  This  Committee  consisted  of: 


t 


Mr. 

C. 

F. 

Corns 

Mr . 

R. 

L. 

Delyea 

Mr. 

R. 

F. 

Malm,  OCE  Coordinator 

Mr . 

L. 

G. 

Guthrie 

Mr. 

D. 

B. 

Baldwin 

Mr . 

R. 

A. 

McMurrer 

The  New  Orleans  District  did  a remarkable  job  in  playing  hosts  to 
the  Conference. 

There  were  13  Division  speakers,  25  moderators,  2 invited  speakers, 
b technical  speakers,  and  10  session  chairman,  who  shared  the  technical 
load  of  the  Conference.  Also,  8 computer  vendors  showed  their  ware  to 
the  participants. 

The  editor  would  like  to  thank  all  the  individuals  who  served 
on  the  committees  and  the  speakers  and  the  moderators  for  sharing 
their  time  and  thoughts.  Without  them  the  Conference  would  not  have 
been  the  success  it  was.  Mr.  Donald  Dressier,  LMVD,  and  Mr.  William 
Price,  WES,  are  specially  thanked  for  their  technical  guidance  and 
assistance . 

This  report  was  edited  by  Dr.  N.  Radhakrishnan,  Research  Civil  En- 
gineer, Computer  Analysis  Branch  (CAB)  and  Special  Technical  Assistant, 
ADP  Center,  under  the  direct  supervision  of  Mr.  J.  B.  Cheek,  Jr.,  Chief, 
CAB,  and  under  the- general  supervision  of  Mr.  D.  L.  Neumann,  Chief, 

ADP  Center. 

The  Director  of  WES  during  the  Conference  and  the  preparation  of 
this  report  was  COL  G.  H.  Hilt,  CE.  Mr.  F.  R.  Brown  was  Technical 
Director . 
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SINGLE-CELL  CONDUITS 


by 

George  W.  Henson* 

Scope 

This  paper  discusses  the  computer-aided  design  of  cast-in-place 
single-barrel  conduits  through  embankment-type  dams  and  levees.  Em- 
phasis is  on  circular  and  oblong  shapes  under  high  fill.  General  pur- 
pose programs  which  could  be  used  will  not  be  dicussed. 

Design  Philosophy 

The  criteria  for  design  of  conduits  used  in  the  Corps  of  Engineers 
(CE)  are  defined  in  EM  1110-2-2902  of  March  3,  1969,  Conduits,  Culverts 
ana  Pipes.  The  EM  recommends  that  circular  shapes  (see  Figure  l)  be 
analyzed  as  a ring  of  uniform  thickness,  since  the  variations  in  thick- 
ness in  the  lower  half  of  the  conduit  may  be  disregarded  in  the  de- 
sign without  appreciable  error.  The  following  is  a summary  of  this  EM. 
Loads 

The  EM  describes  when  loads  are  to  be  used  and  their  magnitude. 
These  loads  are  (a)  groundwater  and  surcharge  water,  (b)  internal  water 
pressure,  (c)  concentrated  live  loads,  and  (d)  backfill. 

Loading  conditions  due  to  backfill 

The  backfill  loads  are  described  for  the  following  construction 
conditions  (see  Figure  2).  Condition  I applies  to  conduits  completely 
buried  in  a ditch  without  superimposed  fill  above  the  top  of  the 
ditch.  Condition  II  applies  to  conduits  completely  buried  in  a ditch 
with  superimposed  fill  above  the  top  of  the  ditch.  Condition  III 
applies  to  conduits  that  project  above  an  embankment  subgrade.  Con- 
dition IV  (not  illustrated)  applies  to  special  conditions  encountered 


* Structural  Engineer,  Tulsa  District. 
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CONDITION  I 


CONDITION  IE 


Figure  2.  Loading  condition  diagrams  (from  EM  1110-2-2902, 
Appendix  III,  3 Mar  1969) 


T 


to  warrant  deviation  from  the  above  loading  conditions. 

The  loads  given  in  Conditions  I and  II  are  similar  to  those  used 
in  industry  for  the  types  of  construction  described.  Conditions  III 
and  IV  are  more  likely  to  describe  construction  conditions  used  in  CE 
outlet  works  structures. 

Design  stresses 

The  allowable  stresses  are  the  same  as  those  used  in  the  ACI  code 
except  for  shear. 

f c ' = 1,000  psi  or  5*000  psi  (except  for  small  conduits  under  low 

fill) 

fc  = 0.15  fc'  (max  fiber  stress) 

fs  = 20,000  psi  allowable  tension  on  intermediate  grade  steel 

Cast-in-place  conduits  - shear.  Shearing  or  diagonal  tension 
stresses  are  not  considered  to  be  as  critical  as  they  are  in  beams  and 
slabs.  The  compression  in  the  ring  at  points  of  maximum  shear  acts  to 
decrease  the  principal  diagonal  tension  stresses,  and  the  diagonal 
tension  requirements  are  less  rigid  than  the  ACI  code.  Shear  and  diag- 
onal tension  requirements  have  been  met  when  the  principal  diagonal 
tension,  at  points  of  maximum  shear,  does  not  exceed  2/fc''  . This 
may  be  determined  by  the  formula: 


where  f^  is  the  diagonal  tension  and  v is  the  average  shearing 
stress. 

Rectangular  box  culverts  - shear.  Shearing  stresses  and  safety 
factors  in  shear  should  be  determined  in  accordance  with  the  University 
of  Illinois  "Development  of  Design  Criteria  for  Reinforced  Concrete  Box 
Culverts,"  Part  II,  "Recommendations  for  Design,"  Structural  Research 
Series  No.  l6U. 

Design  procedure 

General . Several  methods  of  determining  the  maximum  moments, 
thrusts,  and  shears  for  conduits  with  curved  shells  are  given  in  the  EM. 
Both  steel  and  concrete  stresses  at  the  inside  face  of  conduits  with 
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curved  shells  should  be  corrected  for  curved  beam  effects. 

Calculations ■ Cast- in-place  conduits  must  be  designed  for  bend- 
ing moment  combined  with  thrust.  The  EM  recommends  designing  by  the 
transformed-section  method  and  then  investigating  the  critical  sections 
by  ultimate  strength  design  methods.  The  minimum  load  factor  should  be 
1.8,  and  the  capacity  reduction  factor  <)>  should  be  0.7* 

Reinforcement . The  minimum  steel  requirement  of  0.002  bt  (depth 
x thickness)  in  each  face  and  a minimum  cover  of  U inches  are  the  major 
deviations  from  the  ACI  code. 

Computer  Design/Analysis 

The  program  should  be  so  well  written  and  documented  that  a 
structural  engineer  with  little  computer  experience  can  tell  what  the 
program  does  and  how  to  correctly  use  it.  Some  of  the  more  pertinent 
stipulations  are  as  follows: 

Documentation 

The  write-up  should  completely  describe  the  problem  solved  by  the 
program  and  the  restrictions  or  limitations  which  should  be  observed. 

It  should  have  a series  of  statements  describing  the  solution  techniques, 
input  requirements,  output  definitions,  and  any  additional  information 
that  would  be  of  assistance  in  understanding  the  manner  of  solutions  and 
answers  supplied  by  the  program.  An  example  case  with  a detailed  de- 
scription of  all  input  data  and  how  it  is  to  be  prepared  should  be 
included . 

Input 

Problem  identification  is  information  read  and  listed  at  the  begin- 
ning of  the  printout  to  identify  the  problem,  the  user,  or  anything  else 
pertaining  the  run.  Geometry  is  input  describing  the  size  and  shape  of 
the  structure.  The  program  should  generate  most  of  the  dimensions  for 
the  standard  shapes  shown  in  Figure  2 of  EM  1110-2-2902.  When  conduits 
of  constant  thickness  are  being  designed,  the  program  should  also  in- 
crease the  wall  thickness  of  these  shapes  if  they  are  inadequate  in 
shear  or  moment.  Loads  include  the  magnitude  and  location  of  loads 
that  will  be  placed  on  the  structure.  It  should  be  possible  to  check 
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several  loading  cases  on  one  run  without  reentering  the  geometry  and 
design  data.  Design  criteria  include  the  input  of  allowable  stress  or 
safety  factors,  strength  of  concrete  and  steel,  and  minimum  and/or 
maximum  reinforcement  to  be  used  in  the  design.  The  format  statements 
should  be  written  so  that  blanks  or  zeros  cause  the  program  to  insert 
the  most  commonly  used  values  as  inpuc. 

Analys i s 

The  analysis  and  design  should  be  in  accordance  with  the  criteria 
set  forth  in  EM  1110-2-2902.  If  any  other  information  such  as  the  per- 
missible shear  stress,  according  to  the  ACI  code,  is  computed,  it  should 
be  listed  as  information  only  and  not  as  part  of  the  design.  The  program 
should  have  the  option  of  a complete  design  or  of  analyzing  a section 
with  preset  wall  thickness  and  reinforcement. 

Output 

All  output  data  should  have  a neat  and  easily  read  format.  All 
data  should  be  labeled;  and  where  appropriate,  the  units  should  be  in- 
dicated. The  following  items  should  be  included  in  the  output. 

ei.  The  input,  printed  and  labeled. 

Id.  Location  on  the  structure  where  the  moments  and  shears  are 
being  computed. 

c_.  Concrete  thickness. 

d.  Thrust. 

£.  Shear . 

f\  Moment . 

g.  Reinforcing  ratio  or  steel  area. 

h.  Steel  stress. 

_i.  Concrete  stress. 

j_.  Diagonal  tension. 

It.  Points  of  contraflexure  for  box  culverts. 

1_.  Curved  beam  correction  factor. 

m.  Factor  of  safety  in  ultimate  strength  design. 

n.  A printed  summary  (if  more  than  one  case  is  computed  in  a 
single  run)  giving  the  maximum  stresses  at  each  voussoir. 
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Programs  Available 

The  following  programs  have  been  written  by  CE  personnel. 

713C2390  - Omaha  District 

The  Conduit  Design  Program  was  written  in  FORTRAN  IV  with  a re- 
quired machine  storage  capacity  of  approximately  11. 5K.  The  purpose 
of  the  program  is  to  design  single-barrel  conduits.  It  is  also  capable 
of  designing  symmetrical  double-barrel  conduits  with  symmetrical  load- 
ings, but  this  mode  is  not  completely  debugged.  The  program  computes 
geometry,  required  areas  of  reinforcing  steel,  and  stresses.  Two  basic 
modes  of  the  program,  circular  and  rectangular,  are  to  be  considered. 

The  circular  mode  is  the  most  versatile  as  it  encompasses  horseshoe 
sections  and  oblong  sections.  The  rectangular  mode  of  the  program  con- 
siders only  rectangular  openings  with  rectangular  walls.  The  two-barrel 
mode  will  not  be  discussed  except  where  input  data  apply. 

Analysis . The  program  has  three  segments.  The  first  segment  is 
a geometry  phase  that  provides  the  coordinate  description  of  the  ele- 
ments (voussoirs)  for  which  the  analysis  is  made.  The  program  has  the 
capability  of  some  variation  in  geometrical  shape  and  number  of  vous- 
soirs to  be  analyzed.  The  second  segment  is  the  arch  analysis  based 
on  the  examples  set  forth  in  the  PCA  information  sheet — Analysis  of 
Arches,  Rigid  Frames,  and  Sewer  Sections.  This  segment  provides  the 
moment,  thrust,  and  shear  at  the  centerline  of  each  voussoir.  The 
third  segment  is  the  stress  analysis  phase  that  is  based  on  the  design 
criteria  set  forth  in  EM  1110-2-2902. 

Documentation.  Omaha  District  is  revising  the  write-up  to  a more 
usable  form. 

Input.  There  is  a different  number  of  input  cards  for  each  mode 
of  the  program.  The  first  card  contains  fixed  point  variables  that 
control  the  mode  and  the  output.  Sample  data  coding  sheets  and  sketches 
are  used  to  clarify  formats  and  geometry.  Card  No.  1 will  be  defined; 
then  the  input  and  the  output  for  each  mode  will  be  described. 

CARD  NO.  1.  The  first  input  card  contains  fixed  point  variables 
ISENS1 , ISENS2 , ISENS3,  and  ISENSI4. 
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Card 

Column 


Description  Forms 

ISENS1  controls  intermediate  printout  for  debugging  I101 
the  program;  -1  and  0 = YES,  and  +1  = NO;  i.e., 
suppresses  intermediate  printout. 

ISENS2  applies  only  to  the  two-barrel  mode;  +1  = I10j 

straight  leg  and  0 = arch  on  leg.  Use  +1  for  all 
other  modes. 

ISENS3  produces  punch  card  output  for  moment,  I104 

thrust,  and  shear;  -1  and  0 = YES,  and  +1  = NO. 

This  permits  use  of  a second  program  for  design 
with  input  cards  already  punched. 

ISENSU  controls  the  mode.  -1  = rectangular  mode  I10j 

-1  geometry  data  card 

0 = circular  mode  (includes 
horseshoe  and  oblong  sections) 
-4  geometry  data  cards 


1 = two-barrel  mode 

-2  geometry  data  cards 


CIRCULAR  MODE 
Input : 

CARD  NO.  1. 


As  described  above. 


CARDS  NO.  2 thru  NO.  5.  Geometry  cards.  All  dimensions 
are  in  feet.  Note:  All  horizontal  or  X dimensions  are 

negative  (-)  for  input  cards  (except  when  they  are  zero). 
Except  for  three  values  on  Card  No.  2,  all  the  data  are 
defined  in  Figure  3 of  the  user's  manual.  The  three 
variables  not  defined  are  THET,  MID,  and  KOTAL. 


11-20 

THET  = One-half  of  the  incremental  angle  of 
voussoir**  (floating  point  variable).  THET 
in  Figure  3. 

each 

is  noted 

F10.0 

21-30 

MID  = Number  of  voussoirs**  plus  one  (fixed 
variable).  The  maximum  number  of  voussoirs 

point 
is  34 . 

110* 

31-40 

KOTAL  = Twice  the  number  of  voussoirs  plus  one 
(fixed  point  variable). 

110* 

CARD  NO.  6.  Title  card.  Any  alphanumeric  description 
of  a project  may  be  written  in  columns  1 thru  80. 

CARD  NO.  7-  As  described  below. 


* These  values  are  to  be  right  justified. 

**  The  circular  mode  divides  the  conduit  radially  by  incremental 
angles  into  elements  or  voussoirs. 
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Card 

Column 


11-20 


21-30 

31-40 

41-50 


11-20 

21-30 


31-40 

41-50 


51-55 


56-60 


61-70 

71-75 


Description  Format 

NCARDS  = Number  of  loading  cases  for  a given  15* 

geometry.  Cards  8 and  9 are  a loading  case,  and 
any  number  of  sets  of  Cards  8 and  9 may  be  included 
in  the  data  deck. 

LIMIT  = Twice  the  number  of  voussoirs  plus  one  (same  15* 
as  K0TAL). 

BIGK  = 1/2  fckj  as  selected  from  concrete  working  F10.3 
stress  design  handbooks. 

STDEST  = Allowable  steel  stress  in  psi.  F10.3 

DESJ  = j for  working  stress  design.  F10.3 

COVER  = Minimum  concrete  cover  on  reinforcing  steel.  F10.3 
NO.  8.  As  described  below. 

UVINT  = Loading  uniform  vertical  intensity  psf.  F10.0 

UHINT  = Loading  uniform  horizontal  intensity  psf.  F10.0 

THINT  = Loading  triangular  horizontal  intensity  F10.0 

psf  (conduit  height**  wt/cu  ft). 

CVINT  = Loading  concentrated  vertical  at  loc  pounds.  F10.0 


11-20 


21-30 

31-40 


CVINT1  = Loading  concentrated  vertical  at  loci 
pounds. 

LOC  = Selected  location  for  concentrated  load 
(must  be  even  number). 

LOCI  = Selected  location  for  concentrated  load 
(must  be  even  number). 

WTC0NC  = Weight  of  concrete  pcf. 

SYMM  = Blank  for  fixed  end  (two-barrel  mode); 

1 for  symmetrical  single  opening. 

NO.  9.  As  described  below. 

AASUBS  = Initial  area  of  tension  steel  assumed. 

APSUBS  = Initial  area  of  steel  in  compression 
face  assumed. 

FPSUBS  = Concrete  28-day  strength  in  psi. 

XN  = N ratio  of  modulus  of  elasticity  of  steel 
and  concrete. 


F10.3 


F10.3 

F10.3 

F10.3 


F10.0 

F10.0 

F10.0 


F10.0 


* These  values  are  to  be  right  justified. 

**  The  circular  mode  divides  the  conduit  radially  by  incremental 
angles  into  elements  or  voussoirs. 
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Card 

Column 

Description 

Format 

Ul-50 

ASTSTR 

= Allowable  steel  stress  in  psi. 

F10.0 

51-60 

AFSUBC 

= Allowable  concrete  flexure  stress  in  psi. 

F10.0 

61-70 

R = Interior  radius  used  in  computing  curved  beam 
factor. 

F10.0 

71-80 

XINC  = 

Increment  of  increase  in  steel  area. 

F10.0 

Output : 

THE  FIRST  PAGE  of  output  gives  interior  (XI,  YI ) and  exterior  (XO, 
YO)  rectangular  coordinates  along  a radial  line  at  an  angle  ZETA.  ZETA 
is  measured  from  the  Y axis  clockwise.  Note  that  the  origin  of  ZETA  is 
not  the  same  as  the  rectangular  coordinate  axis.  The  voussoirs  are  cre- 
ated by  an  incremental  angle  radially  from  the  center  of  the  opening. 

THE  SECOND  PAGE  is  the  input  geometry. 

THE  THIRD  PAGE  gives  interior  (XI,  Yl),  exterior  (XO,  YO),  and 
centerline  (XC,  YC)  rectangular  coordinates. 

THE  FOURTH  PAGE  of  output  gives  three  values.  The  total  voussoir 
area  is  the  total  cross-sectional  area  of  the  conduit.  VTOTAL  is  one- 
half  the  total  vertical  pressure  (pounds).  UBASEP  is  the  uniform  base 
pressure  in  pounds  per  square  foot. 

THE  FIFTH  PAGE  of  output  prints  the  input  loads,  then  lists  the 
moment  thrust  and  shear  at  the  centerline  of  each  section.  The  actual 
thickness  of  each  voussoir  is  the  next  column.  The  following  columns 
are  the  eccentricity  of  the  axial  load,  the  maximum  and  minimum  con- 
crete stresses,  and  the  diagonal  tension. 

THE  SIXTH  PAGE  of  output  shows  the  required  d of  the  concrete  for 
balance  design  (WSD).  The  next  column  gives  the  actual  d provided. 
Following  these  columns  are  the  Fixed  End  Movement,  Thrust  at  Fixed  End, 
and  Shear  at  Fixed  End.  These  are  useful  only  in  the  two-barrel  mode. 

THE  SEVENTH  PAGE  of  output  prints  on  the  first  line  input  material 
specifications.  Below  follows  the  concrete  and  reinforcing  steel 
stresses;  also  the  required  areas  of  reinforcing  steel.  The  input 
minimum  steel  area  was  that  required  for  temperature  and  shrinkage 
stresses.  Negative  signs  in  the  column  TENSION,  STEEL  (PSI)  indicates 
that  iteration  for  kd  (WSD)  did  not  converge  and  erroneous  value  was 

Ik 


used;  therefore  all  values  for  that  voussior  are  not  correct.  This  usu- 
ally occurs  when  tension  in  the  tension  face  of  voussior  is  very  small. 
RECTANGULAR  MODE 


r 


Card 

Input 

Column 

CARD 

CARD 

Description 

NO.  1 is  described  earlier. 

NO.  2 is  the  geometry  card  (Figure  4). 

Format 

1-6 

AA 

F 6.0 

7-12 

AB 

F6.0 

8-l8 

AC 

F6. 0 

19-2U 

AD 

f6.o 

25-30 

AF 

f6.o 

31-36 

AG 

f6.o 

37-^0 

DIVAA 

FU.O 

Ui-UU 

DIVAB 

F^  .0 

l»5-l»8 

DIVAC 

Fk.O 

L9-52 

DIVAD  = number  of  elements  in  segment  AD 

FU.O 

53-56 

DIVAF 

FU.O 

57-60 

DIVAG 

FU.O 

61-65 

CARDS 

NO  = 2*  (DIVAF  + 2*  DIVAG  = 2*  DIVAD)  + 1 
3 thru  5 are  omitted. 

15 

CARDS  6 thru  9 are  the  same  as  described  for  the  circular  mode. 
On  Card  No.  7,  LIMIT  is  now  equal  to  NO. 


Output : 

The  first  page  of  output  gives  interior  (XI,  Yl)  and  exterior 
(XO,  YO)  rectangular  coordinates  of  each  voussoir  and  the  centerline  of 
each  voussoir.  There  is  no  angle  ZETA. 

The  rest  of  the  output  is  described  by  the  write-up  for  the  cir- 
cular mode.  The  only  difference  in  results  is  that  the  voussoirs  are 
irregularly  defined  rather  than  created  by  an  incremental  angle  radially 
from  the  center  of  the  opening. 

713-G1-M0090  - Southwestern  Division 

This  program  computes  moments,  thrusts,  shears,  and  factors  of 


L*  H> I - - I ■ Ml— II  II  M 


AF  & DIVAF 


safety  in  shear  and  combined  axial  load  and  flexure  for  conduits  under 
high  fill.  It  designs  or  analyzes  circular,  oblong,  and  horseshoe  shapes 
of  constant  thickness.  Conduits  with  a variable  cross  section  may  be 
analyzed  if  the  voussoir  coordinates  describing  the  shape  are  input. 

Analysis.  For  the  standard  conduits  of  constant  thickness,  the 
designer  has  the  option  of  analyzing  a specific  situation  or  of  using 
the  program  to  automatically  increment  the  concrete  thickness  and 
percentage  of  reinforcement  to  obtain  the  specified  minimum  factors  of 
safety.  The  program  consists  of  the  following  phases: 

a.  Geometry  phase.  The  shapes  of  the  conduits  of  constant  thick- 
ness are  defined  in  the  input  by  inside  radius,  crown  thick- 
ness, and  by  the  tangent  length  in  the  case  of  an  oblong  con- 
duit. In  the  geometry  phase  the  voussoir  coordinates  are 
computed  for  these  shapes.  Since  9 degrees  was  selected  as 
the  angle  increment,  a circular  conduit  (half  section)  will 
have  20  voussoirs.  A horseshoe  shape  has  21  voussoirs.  An 
oblong  shape  has  20  plus  the  number  of  divisions  input  for  the 
tangent  length.  A variable  cross  section  may  have  any  number 
up  to  30,  but  the  coordinates  will  be  part  of  the  input. 

b.  Load  determination  phase.  Loading  variables  include  the  ver- 
tical and  horizontal  soil  pressure  coefficients,  conduit  in- 
vert elevation,  water  table  elevation,  soil  elevation,  rock 
line  elevation,  and  the  unit  weight  of  moist  and  saturated 
soil.  If  live  load  is  a consideration,  it  may  be  expressed 

as  a surcharge  by  modifying  the  fill  elevation  and/or  the  soil 
pressure  coefficients.  The  effect  of  placing  the  rock  line 
above  the  invert  of  the  conduit  is  to  eliminate  the  horizontal 
soil  loads  below  the  rock  line  elevation.  Soil  and  water 
pressure  are  computed  at  each  voussoir  level,  and  the  effects 
of  conduit  dead  load  are  included.  Uplift  and  base  pressure 
are  computed,  and  the  forces  on  the  voussoirs  are  subscripted 
and  stored  for  use  in  calculations  of  moment,  thrust,  and 
shear  in  the  next  phase. 

c_.  Moment,  thrust,  and  shear  calculation  phase.  The  program  uses 
the  elastic  center  method  of  analysis  for  computing  moments, 
thrusts,  and  shears.  For  a detailed  discussion,  the  reader 
is  referred  to  PCA  pamphlet  ST  53,  Analysis  of  Arches,  Rigid 
Frames,  and  Sewer  Sections.  The  moment  correction  discussed 
in  the  PCA  pamphlet  was  used  in  order  to  yield  more  realistic 
values  for  moments  throughout  the  conduit. 

d_.  Ultimate  strength  phase.  This  phase  of  the  program  computes 
factors  of  safety  in  shear,  and  combined  axial  load  and  flex- 
ure, based  on  the  principles  of  Ultimate  Strength  Design  in 
ACI  318-63.  The  depth  from  the  face  of  the  concrete  to  the 
centroid  of  the  tensile  steel  is  assumed  to  be  5 inches  at 
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any  point.  If  a design  (rather  than  an  analysis)  is  being 
performed,  the  factor  of  safety  in  shear  is  compared  with  the 
factor  of  safety  required  by  the  input;  if  the  computed  value 
is  less  than  the  required  value,  a new  concrete  thickness  is 
computed  and  the  program  returns  to  the  geometry  phase.  This 
computation  is  made  only  for  the  conduits  of  constant  thick- 
ness, and  the  minimum  thickness  increase  is  set  at  3 inches. 

It  should  be  noted  that  the  shear  is  checked  by  the  formulas 
given  in  the  ACI  code  and  not  those  given  in  the  EM.  Under 
ordinary  loading,  this  will  give  thicker  sections  than  would 
be  required  by  the  EM.  The  factor  of  safety  in  flexure  is 
compared  as  it  was  in  shear;  and  if  inadequate,  the  minimum 
percentage  of  reinforcement  input  will  be  increased  by  0.005 
and  the  flexural  computations  repeated.  If  the  flexural 
factor  of  safety  is  adequate  or  if  the  maximum  percentage  of 
reinforcement  input  is  reached,  the  output  will  be  printed. 
While  the  program  increments  the  wall  thickness  until  a sec- 
tion adequate  in  shear  is  obtained,  the  upper  limit  of  p 
may  cause  the  computations  to  stop  at  a point  where  the  flex- 
ural factor  of  safety  is  still  inadequate. 

In  the  analysis  mode,  the  percentage  of  reinforcement  and 
the  wall  thickness  will  remain  constant.  If  more  than  one 
configuration  of  reinforcement  is  desirable,  an  analysis 
should  be  made,  using  each  different  percentage.  In  this 
case,  the  designer  should  realize  that  the  factors  of  safety 
are  valid  only  at  the  voussoirs  where  the  indicated  steel 
ratios  actually  occur. 

Documentation.  The  write-up  is  very  good.  It  describes  the  solu- 
tion routine  and  gives  references.  The  input  requirements  and  the  out- 
put descriptions  are  very  clearly  listed.  Several  example  runs  are 
attached.  Sense-switch  settings  are  a little  vague. 

Input.  A sense  switch  is  required  to  control  selection  between 
design  and  analysis.  It  would  be  better  if  a control  card  had  been  used 
for  this  function. 

Conduits  of  constant  thickness.  Five  cards  are  required  for  a run. 
Any  change  in  loading  or  design  criteria  requires  an  additional  run. 
Twenty-four  items  of  input  data  are  required  for  the  smallest  possible 
input.  The  following  is  a summary  of  the  cards. 

Card  Format  Probable 

Column  Description  Specification  Magnitude 

1-5  CARD  WO.  1 - Header  information.  The  first  TO  columns  of  this 

card  are  used  for  problem  identification.  Information  con- 
tained in  the  card  will  be  read  and  listed  at  the  beginning 


* 


Card 

Column 


1-5 

6-10 

11-15 

16-20 

21-30 

31-10 

Ul-50 

51-60 


Format  Probable 

Description  Specification  Magnitude 

of  the  printout.  The  last  10  columns,  if  used,  will  be 
ignored  by  the  computer.  In  the  IBM  1620  version  of  the 
program,  this  record  is  entered  under  a Hollerith  specifica- 
tion. In  the  GE  225  version,  it  is  entered  under  series 
"A"  specifications  (11A6,  Al). 


CARD  W0.  2 - Geometry  data.  The  second  record  of  input  data 
contains  the  following  eight  values. 


Number  of  voussoirs  (NY)* 

15 

Number  of  divisions  of  tangent 
length  in  oblong  section  (NDIV)** 

15 

Control  parameter  denoting  type  of 
standard  section  (M)t 

15 

Control  parameter  denoting  a 
variable  section  (MM)tt 

15 

Thickness  of  section  at  the  crown, 
in  feet  (TC) 

F10.0 

Inside  radius  of  standard  section, 
in  feet  (RI)* 

F10.0 

Tangent  length  of  oblong  section, 
in  feet  (TANL)** 

F10.0 

Invert  elevation,  in  feet  (ELINV), 
positive  only 

F10.0 

XX 

X 

( + )X 

(-)X 

{+)X 

(-r 

X.XXX 
X.XXX 
X.XXX 
XXX. XX 


* For  circular  conduit,  NY  = 20  . For  horseshoe  conduit,  NY  = 21  . 
For  oblong  conduit,  NV  = 20  + NDIV  . For  nonstandard  section, 

NV  = any  number  not  exceeding  30  . 

**  For  any  type  of  section  other  than  oblong,  NDIV  = 0 . May  be  left 
blank  if  data  is  to  be  processed  by  GE  225. 

t For  circular  conduit,  M = -1  . For  horseshoe  conduit,  M = +1  . 
For  any  other  type  of  section,  M = 0 . May  be  left  blank  if  data 
is  to  be  processed  by  GE  225. 

++  For  variable  section,  MM  = any  number  other  than  zero.  For 

standard  section,  MM  = 0 . May  be  left  blank  if  data  is  to  be  pro- 
cessed by  GE  225. 

♦ For  a nonstandard  section,  RI  = 0 . May  be  left  blank  if  data 
is  to  be  processed  by  GE  225. 

**  For  anything  other  than  an  oblong  section,  TANL  = 0 . May  be 
left  blank  if  data  is  to  be  processed  by  GE  225. 
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L 


Card 

Column 


Format 

Specification 


Probable 

Magnitude 


Description 


CARD  NO.  3 - Soils  data.  The  third  record 
lowing  soils  data  elements. 

consists 

of  the  f ol— 

1-10 

Coefficient  of  vertical  earth 
pressure  (VK) 

F10.0 

X.XX 

11-20 

Coefficient  of  horizontal  earth 
pressure  (HK) 

F10.0 

.XX 

21-30 

Soil  elevation,  in  feet  (SELEV)* 

F10.0 

XXX. XX 

31-^0 

Water  elevation,  in  feet  (WELEV)* 

F10.0 

XXX.  XX 

Ul-50 

Rock  line  elevation,  in  feet  (RELEV)* 

F10.0 

XXX. XX 

51-60 

Unit  weight  of  moist  soil,  in  kef 

F10.0 

.XX 

61-70 

Buoyant  unit  weight  of  saturated  soil, 
in  kef  (WES) 

F10.0 

.XX 

CARD  NO.  4 _ Design  criteria.  The  fourth  record  contains  the 
following  seven  items. 

1-10 

Unit  weight  of  concrete  in  kef 

(wc) 

F10.0 

.XX 

11-20 

Minimum  percentage  of  reinforcement 
in  tension  face  (PMIN) 

F10.0 

.XX 

21-30 

Maximum  percentage  of  reinforcement 
in  tension  face  (PMAX) 

F10.0 

.XX 

31-UO 

28-day  compressive  strength  of 
concrete,  in  psi  (FC)** 

F10.0 

xxxx. 

1*1-50 

Yield  strength  of  reinforcing 
steel,  in  psi  (FY) 

F10.0 

xxxxx. 

51-60 

Required  factor  of  safety  in 
shear  (FSSR)t 

F10.0 

X.XX 

61-70 

Required  factor  of  safety  in 
flexure  (FSFR)t 

F10.0 

X.XX 

CARD  NO.  5 - Design  criteria.  The  fifth  record  of 
contains  only  two  data  elements,  as  follows. 

input 

1-10 

Ultimate  strength  capacity  reduction 
factor  for  shear  (PHIS) ft 

F10.0 

.XX 

* All  elevations  must  be  entered  as  positive  values. 


**  Compressive  strength  of  concrete  (FC)  = i+000  psi  minimum. 

t Usually  1.8. 

++  Usually  0.85. 
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Card 

Column 


Format 
>ecif  ication 


Probabl e 
Magnitude 


Description  Sp 

11-20  Ultimate  strength  capacity  reduction  F10.0  .XX 

factor  for  combined  axial  load  and 
flexure  (PHIF)* 

RKMAIiiLHU  CARD.'  - Voussoir  definition  data.  The  remaining  cards 
of  the  input  deck  are  used  to  define  the  shape  of  a variable- 
section  conduit  and  are  therefore  not  required  for  the  standard 
conduits  of  constant  thickness.  There  should  be  one  card  for 
each  voussoir  of  the  variable  section,  and  each  card  contains 
four  data  elements  as  follows. 


1-10 

X-coordinate 
in  feet  (X^ ) 

of 

voussoir 

extrados , 

F10.0 

X.XXX 

11-20 

Y-coordinate 
in  feet  (Y.) 

of 

voussoir 

extrados , 

F10.0 

X.XXX 

21-30 

X-coordinates  of  voussoir 
in  feet  ( Z . ) 

intrados , 

F10.0 

X.XXX 

31-40 

Y-coordinate 
in  feet  (Vf  ) 

of 

voussoir 

intrados , 

F10.0 

X. 

Output . All  output  from  this  program  is  produced  on  the  on-line 


printer  (or  on  the  console  typewriter,  if  being  run  on  the  IBM  1620). 
For  purposes  of  visual  verification,  the  basic  input  data  is  listed 
along  with  the  computed  output  values.  The  input  items  are  printed  in 
the  same  sequence  and  format  (i.e.,  one  line  per  card)  in  which  they 
are  entered,  and  are  identified  by  their  mnemonic  labels.  The  voussoir 
coordinates — whether  computed  or  entered  as  input  data — are  listed  fol- 
lowing printout  of  the  basic  input  data.  Output  of  the  computed  values 
(excluding  voussoir  coordinates)  is  as  follows: 

Voussoir  number  (beginning  at  zero) 

Voussoir  thickness  (inches) 

Percentage  of  reinforcement 

Thrust  (axial  load)  acting  at  center  of  voussoir  (pounds) 

Shear  acting  at  center  of  voussoir  (pounds) 

Final  moment  due  to  all  loading  (inch-pounds) 


Note:  All  voussoir  coordinates  have  their  origin  at  crown  extrados  and 

must  be  entered  as  positive  values. 

* Ususally  0.70. 
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Factor  of  safety  in  shear 
Factor  of  safety  in  flexure 
i -S-K2-H223  - Louisville  District 

The  program  was  written  to  design  a variable  seciton  conduit  with 
water  and/or  earth  loading.  This  is  a revised  version  of  Tulsa  District 
program  13-G1-G523  written  in  GE  WIZ  for  a GE  225  computer.  The  Load 
Determination  Phase  and  the  Moment,  Thrust,  and  Shear  Phase  of  the 
Southwestern  Division  program  were  also  taken  from  the  Tulsa  District 


program. 


Analysis.  The  cross  section  of  the  conduit  is  divided  into  a 


maximum  of  30  wedge-shaped  pieces  by  inputing  the  coordinates  of  the 
corners  of  the  voussoirs.  All  moments,  shears,  and  thrusts  are  com- 
puted by  the  elastic  center  method  on  a line  bisecting  the  outside  and 
inside  edge  of  each  voussoir.  The  external  loads  are  applied  on  free- 
bodies  at  the  outside  face  of  the  conduit.  Moments  are  computed  to  the 
center  of  horizontal  and  vertical  loads.  Moments  are  then  summed  about 
the  center  of  the  voussoirs.  This  is  opposed  to  the  PCA  type  of  anal- 
ysis where  all  load  is  assumed  to  be  concentrated  at  the  center  of  each 
voussoir.  The  procedure  used  in  the  program  produces  lower  moments 
than  other  methods,  but  the  results  may  be  more  accurate. 

Loading  may  be  varied  by  adjusting  the  height  of  water,  saturated 
soil  or  moist  soil;  by  adjusting  the  weight  of  saturated  or  moist  soil; 
or  by  adjusting  the  coefficients  of  horizontal  and  vertical  earth  pres- 
sure. The  horizontal  external  soil  load  is  eliminated  below  the  rock 
line  elevation.  This  program  corrects  stresses  at  the  interior  face 
for  curved  beam  effects. 

Documentation.  The  write-up  is  the  same  as  the  Tulsa  District 
program.  It  should  be  revised  to  include  a description  of  the  curved 
beam  correction  factor  and  how  it  is  to  be  entered  as  input. 

Input . The  input  described  below  is  good  for  one  loading  case. 
Additional  cases  would  require  another  complete  set  of  input  data. 


Card 

Column 


l Description 

CARD  NO.  1.  Project  identification  card. 


r 


Card.  Format 

Column  Description  Specification 


CARD  NO.  2. 


2-8 

(CC)  Distance  to  tensile  steel  (inches) 

F7.2 

9-16 

(HD)  Height  of  water  above  (+)  or  be- 
low (-)  extrados  at  top  of  transition  (feet) 

F8.2 

17-18 

(HES)  Height  of  submerged  soil  above  (+)  or 
below  (-)  extrados  at  top  of  transition  (feet) 

F8.2 

25-32 

(WES)  Unit  weight  of  submerged  soil  (kef) 

F8.4 

33-1*1 

(HEM)  Height  of  moist  soil  above  HES  if 
HES  is  ( + ).  Above  or  below  crown  if  HES 
is  (-)  (feet) 

F8.2 

Ul-48 

(WEM)  Unit  weight  of  moist  soil  (kef) 

F8.4 

1*9-56 

(RL)  Rock  line  distance  below  extrados  at  top 
of  transition  (feet) 

F8.2 

57-64 

(KH)  Coefficient  of  horizontal  soil  pressure 

F8.2 

65-72 

(KV)  Coefficient  of  vertical  soil  pressure 

F8.2 

73-80 

(NR)  Modulus  of  elasticity  ratio 

F8.2 

CARD  NO.  3. 

2-8 

FS  - Allowable  steel  working  stress  (ksi) 

F7.3 

9-16 

FC  - Ultimate  concrete  compressive  strength 
(ksi) 

F8.3 

17-24 

WC  - Unit  weight  of  concrete  (kef) 

F8.4 

25-32 

IV  - Number  of  voussoirs  (right  justified) 

18 

33-40 

TC  - Thickness  of  crown  (feet) 

F8.2 

4i-48 

VC  - Index  of  compressive  strength  of  concrete 

F8.4 

49-56 

IS  - Continuous  ends  = 1;  fixed  ends  = 0 
(right  justified) 

18 

CARD  NO.  4.  There  should  be  one  card  for  each  voussoir  of  the 
section,  and  each  card  contains  five  data  elements  as  follows. 

2-8 

I - Voussoir  number  (right  justified) 

17 

9-16 

X(l)  - Abscissa  of  extrados  (feet) 

F8. 2 

17-24 

Y(l)  - Ordinate  of  extrados  (feet) 

F8. 2 

25-32 

Z(l)  - Abscissa  of  intrados  (feet) 

F8.2 

33-40 

W(l)  - Ordinate  of  intrados  (feet) 

F8.2 

LAST  CARD.  Curved  beam  correction  factor. 

2-3 

Number  of  curves 

12 

1 

i 
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Card 

Column  Description 

4-11  Radius  of  curvature  at  top 

12-13  Voussior  number  beginning  curve 

14-15  Voussoir  number  ending  curve 

16-23  Radius  of  curvature  at  bottom 

24-25  Voussoir  number  beginning  curve 

26-27  Voussoir  number  ending  curve 


Format 

Specification 

F8.0 

12 

12 

F8.0 

12 

12 


Output . All  of  the  input  data  are  printed  for  references.  The 
following  computed  values  are  output  for  each  voussoir. 

Voussoir  number  (beginning  at  zero) 

Final  moment  (kip-ft) 

Shear  (kips) 

Thrust  (kips) 

Concrete  stress  (ksi) 

Steel  stress  (ksi) 

Area  of  steel  (square  inches) 

Diagonal  tension  (ksi) 

Curvature 

Curved  beam  correction  factor 

The  weight  in  kips  and  the  volume  in  cubic  yards  of  a 1-foot- 
thick  section  described  by  the  voussoir  points  are  printed  at  the  end 
of  each  run. 

713-G9AU-06q  - Vicksburg  District 

The  program  was  written  to  design  or  analyze  single-cell  box  cul- 
verts. It  determines  the  minimum  thicknesses  of  the  horizontal  and  ver- 
tical members  and  the  area  of  reinforcing  steel  to  provide  for  moment 
and  the  required  factor  of  safety  for  cracking  load  for  shear.  For  a 
maximum  of  30  loading  conditions,  the  program  will  find  the  minimum 
thicknesses  required  to  support  the  load  or  it  will  compute  the  resulting 
stresses,  d's  , steel  reinforcing,  and  factors  of  safety  in  the  mem- 
bers. The  program  does  not  compute  for  haunches.  Backfill  is  con- 
sidered to  be  level  over  the  completed  structure.  Concrete  cover  for 
reinforcement  is  2.5  inches  in  all  cases  except  the  bottom  of  the 
bottom  slab  where  3.5  inches  are  used. 
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Analysis . The  design  is  in  accordance  with  the  criteria  set  forth 
in  EM  1110-2-2902.  Distributed  moments  are  adjusted  to  the  column  face 
and  at  the  center  of  the  span  by  the  method  explained  starting  on  page 
28  of  Continuity  in  Concrete  Building  Frame,  Practical  Analysis  for 
Vertical  Load  and  Wind  Pressure,  fourth  edition,  published  by  the  Port- 
land Cement  Association. 

a.  The  unit  loading  on  the  structure  is  computed  at  the  begin- 
ing  of  computation  of  each  case  and  each  time  there  is  a 
change  in  the  thickness  of  any  of  the  members.  Vertical 
weight  of  the  soil  is  controlled  by  vertical  load  factors  for 
each  case.  A full  culvert  is  used  to  compute  the  forces  for 
internal  load. 

b_.  In  the  design  phase  (no  member  thickness  input),  the  safety 
factor  for  cracking  shear  load  is  used  to  find  the  thickness 
of  the  top  slab.  The  required  safety  factor  (2  to  2.5)  is 
based  on  the  L/d  ratio.  The  sidewall  thickness  is  set  equal 
to  the  top  slab,  and  the  bottom  slab  is  set  1 inch  thicker 
than  the  top  slab.  The  minimum  slab  thickness  is  9 inches. 

c_.  Fixed  end  moments  and  moment  distributions  are  made  between 
the  midpoints  of  the  diagonals  at  the  corners  of  the  struc- 
ture. Distributed  moments  are  adjusted  to  the  column  face. 

d_.  Computations  are  made  to  find  the  unit  loads,  moments,  reac- 
tions, shears,  thrust,  and  required  d's  for  the  structure, 
from  external  forces,  then  internal  forces,  and  finally  hydro- 
static forces.  In  the  design  phase,  the  required  d's  are 
compared  with  the  provided  d's  at  the  end  of  each  force 
computation.  If  the  required  d is  more  than  0.1  inch 
greater  than  the  provided  d , the  member  thickness  is  in- 
creased 1 inch  and  the  computations  are  repeated  beginning 
with  the  external  forces. 

e_.  At  this  point  the  provided  d's  equal  (or  are  less  than  0.1 
inch  smaller)  or  exceed  the  required  d's  , Final  computa- 
tions or  output  are  as  follows: 

(1)  Search  for  maximum  shear  at  the  face  of  each  structural 
member . 

(2)  Maximum  moment,  thrust,  and  Nd"/12  at  each  point  used  to 
compute  the  maximum  steel  moment. 

(3)  Maximum  steel  moment. 

(U)  Required  and  provided  d's  at  each  point. 

(5)  Equivalent  areas  of  steel  for  the  steel  moment  and  the 
thrust . 

(6)  Area  of  steel  reinforcing  required.  This  area  will  not 
be  less  than  the  minimum  factor  times  12d. 
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(7)  Di  stances  between  the  points  of  contraflexure  of  top  and 
bottom  slabs  for  each  maximum  and  external  moment. 

(8)  True  safety  factor  and  required  safety  factor  for  cracking 
load  for  shear  in  the  top  and  bottom  slabs.  The  required 
d does  not  always  control  the  thickness  to  be  used. 

This  is  the  reason  for  computing  the  true  safety  factor 
at  the  end  of  each  computation. 


After  the  starting  thicknesses  are  computed,  if  the  provided  d 
is  1 inch  larger  than  the  largest  required  d of  a member,  the  thick- 
ness is  reduced  1 inch  and  the  computations  for  the  structure  are  re- 
peated. If,  however,  the  provided  d does  not  contain  the  largest  re- 
quired d of  a member  or  does  not  meet  safety  factor  requirements,  a 
flag  is  set  so  that  the  member  cannot  be  reduced  again,  the  thickness 
is  increased  1 inch,  and  the  computations  for  the  structure  are  re- 
peated. In  the  process  of  finding  the  structure  with  the  least  required 
safety  factor  and  the  minimum  required  d , several  safe  combinations 
of  thicknesses  may  be  found.  For  each  safe  combination,  the  cross- 
sectional  area  is  computed.  The  safe  combination  with  the  minimum  end 


area  will  be  used  in  calculating  the  next  case.  These  thicknesses  can- 
not be  reduced  by  subsequent  cases;  however,  they  can  be  increased  if 
the  loading  conditions  require  more  strength  for  moment  or  shear. 

Documentation.  The  write-up.  is  very  good.  It  describes  the  solu- 
tion routine  and  gives  references.  The  input  requirements  and  the  out- 
put descriptions  are  very  clearly  listed.  An  example  problem  is 
attached. 

Input . 

Ei.  USER  ID.  The  first  card  of  the  input  is  to  inform  the  com- 
puter operator  who  requested  the  run  and  where  to  send  the 
output.  Computer  output  is  often  saved  on  disks  and  is  out- 
put by  the  printer  at  a later  time.  Thus,  the  output  should 
contain  the  user  identification  to  ensure  that  it  is  returned 
to  the  sender.  This  card  should  contain  the  requester's  name, 
the  office  where  he  works,  and  the  district  office;  e.g.. 

Bill  Jones,  Minor  Struct  Sec,  VXD. 

Id.  TITLE  CARD.  The  title  card  is  the  second  type  of  input.  It 
contains  80  alphanumeric  characters  to  identify  the  culvert 
being  designed.  The  contents  of  this  card  will  be  printed 
at  the  top  of  each  sheet  of  output. 
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Card 

Column 

6-10 
11-15 
16-20* 
21-25* 
26-30* 
31-35** 
36-Uo 
1*1-1*  5+ 
1*6-50 
51-55++ 
56-79 
80 


COMMENT  CARD.  A comment  card  is  the  third  type  of  input. 

The  first  79  card  columns  contain  the  comment  (in  alpha- 
numerics);  and  card  column  80  must  contain  a minus.  Any 
number  of  comment  cards  may  be  used  between  the  Title  Card 
and  the  Constant  Card.  This  card  is  optional.  It  may  be 
used  to  add  some  special  comment  concerning  the  culvert 
design  for  later  reference.  The  comments  are  printed  under 
the  title  on  the  User  ID  sheet  that  is  printed  for  the  op- 
erator. When  two  or  more  culverts  are  designed,  comments 
are  printed  on  the  sheet  between  the  summary  sheet  of  the 
preceding  culvert  and  the  first  case  of  the  next  culvert. 

CONSTANTS.  The  fourth  type  of  input  card  contains  constants 
that  will  be  used  in  computing  all  of  the  cases.  The  first 
two  entries,  width  and  height,  are  required;  but  the  remainder 
of  the  card  may  be  left  blank.  The  thicknesses  of  the  top, 
side,  and  bottom  are  left  blank  for  design  applications;  but 
they  must  be  entered  when  an  analysis  is  to  be  computed. 
Standard  design  factors  are  entered  by  the  program  when  the 
remainder  of  the  card  is  left  blank. 


Explanation 

Width  of  culvert,  feet. 

Height  of  culvert,  feet. 

Thickness  of  top  slab,  inches. 

Thickness  of  side,  inches. 

Thickness  in  bottom' slab,  inches, 
f^  for  concrete,  psi. 

K factor  for  concrete, 
f for  steel,  psi. 
j factor  for  steel. 

Fraction  of  bd  for  minimum  steel. 

Blank. 

If  iterations  of  the  design  and  backup  computations  are 
desired  in  the  output,  punch  a C in  card  column  80. 


* When  slab  thicknesses  are  entered,  the  design  phase  will  not  take 
place.  Also,  no  summary  sheet  will  appear  at  the  end  of  the  output. 


#* 

Blanks  or 
for  K. 

zero 

cause 

the 

program 

to 

insert 

1*,000 

for 

f ' 

c 

and  32k 

t 

Blanks  or 
0.860  for 

zero 

J. 

cause 

the 

program 

to 

insert 

20,000 

for 

f 

s 

and 

ft 

Blanks  or 

zero 

cause 

the 

program 

to 

insert 

0.0025 

for 

the 

minimum 

steel  factor. 


ji.  CASE  TITLE.  A case  title  card  must  precede  each  case  data 

card.  This  card  uses  the  same  format  as  the  title  card  to  in 
put  up  to  80  alphanumeric  characters.  This  data  is  used  to 
label  the  cases  in  the  output.  It  is  printed  as  the  second 
line  of  input  for  each  case. 

_f.  CASE  DATA.  The  fifth  type  of  input  card  gives  the  input  used 
to  compute  the  loading.  Input  consists  of  groundwater,  as 
many  as  three  soil  conditions,  and  hydrostatic  head. 


Card 

Column 

1-3 

U-5 

6-10 

11-15 

16-20 

21-25 

26-30 

31-35 

36-1+0 

41-1+5 

1+6-50 

51-55 

56-60 

61-65 


Explanation 

The  letters  END  are  to  be  punched  in  the  last  case  only. 
Blank. 

Depth  of  groundwater  in  feet  above  the  top  of  the  opening 
Groundwater  may  extend  above  the  ground  surface. 

Depth  of  soil  in  feet  above  the  top  of  the  opening. 

Weight  of  soil  in  pcf.  Submerged  soils  do  not  include 
the  weight  of  water. 

K factor  for  soil. 

Depth  of  soil  in  feet  above  soil. 

Weight  of  soil,  pcf. 

K factor  for  soil. 

Depth  of  soil  in  feet  above  soil. 

Weight  of  soil,  pcf. 

K factor  for  soil. 

Height  of  head  ic  feet  above  the  top  of  the  opening. 

Each  case  that  contains  head  is  recomputed  as  a separate 
case  with  no  head. 

First  factor  for  computing  the  vertical  weight  of  soil 
on  the  top  slab. 


66-70  Second  factor  for  computing  the  vertical  weight  of  soil 

on  the  top  slab.  This  is  left  blank  except  on  the  rare 
occasion  when  it  is  desirable  to  compute  a case  using  two 
factors.  This  causes  a second  case  computation  to  be 
made. 

Spaces  not  required  for  input  may  be  left  blank.  Each  case  data  card 
must  follow  a case  title  card. 


INPUT  SETUP  Input  for  this  program  must  be  in  the  following 
order. 

(l)  User  ID.  One  time  only. 
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(2)  Title  card.  One  per  culvert. 

(3)  Comment  card.  Optional. 

(k)  Constants.  One  per  culvert. 

(5)  Case  title.  One  per  case. 

(6)  Case  data.  One  per  case. 

The  input  setup  should  include  as  many  pairs  of  case  title  and  case  data 
cards  as  there  are  loading  conditions.  The  last  data  card  must  have 
"END”  in  card  columns  1-3.  The  program  is  compiled  for  a maximum  of 
30  cases.  When  two  or  more  culverts  are  to  be  processed  in  the  same 
run,  the  last  card  of  the  first  culvert  is  followed  by  the  new  title 
card  (b  above).  The  input  setup  after  the  title  card  is  the  same  as  it 
is  for  the  first  culvert. 

Output 

The  information  on  the  User  ID  card  is  the  first  line  of  output 
and  is  followed  by  the  program  number  which  references  the  program  used. 
The  next  item  printed  is  the  title  of  the  culvert.  If  there  are  any 
comment  cards,  they  will  be  listed  after  the  culvert  title.  When  two 
or  more  culverts  are  computed  in  the  same  run,  the  title  of  the  next 
culvert  and  any  comment  cards  to  be  listed  are  printed  on  the  sheet  be- 
tween the  summary  sheet  of  the  preceding  culvert  and  the  first  case  of 
the  next  culvert. 

The  remainder  of  the  output  format  except  the  summary  is  re- 
peated for  each  case.  Each  case  is  output  one  time  for  each  vertical 
load  factor.  If  there  are  10  cases,  two  of  which  have  head,  there  will 
be  12  pages  of  case  output.  Each  of  the  two  having  head  will  be  recom- 
puted for  no  head. 

a.  The  first  line  of  case  output  is  the  culvert  title  and  data 
of  computation.  This  line  is  followed  by  the  case  title. 

_b.  The  first  data  output  are  the  thicknesses  of  the  top,  side, 

and  bottom  members,  input  constants,  and  case  variables.  This 
data  has  headings  for  identification. 

c_.  Output  for  the  external,  internal,  and  hydrostatic  head  load- 
ing uses  the  same  format.  Adjusted  moments  (M),  thrust  (N), 
and  shears  (V)  for  each  loading  are  listed  for  points  1-7. 
These  are  followed  by  maximum  shear,  working  moment,  thrust 
accompanying  working  moment,  equivalent  moment  (ND"/12), 
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steel  moment,  required  d , provided  d , area  of  steel  for 
moment,  equivalent  area  of  steel  for  thrust  (N/f  ),  and  area 
of  reinforcing  steel  required  in  square  inches  for  each  point. 

The  items  above  that  have  separate  calculations  for  inside  and/or  out- 
side faces  have  two  line  entries  with  the  inside  face  on  top. 

d_.  The  safety  factor  and  required  safety  factor  for  the  cracking 
load  for  shear  and  the  distances  between  the  points  of  contra- 
flexure  for  maximum  and  external  loading  are  printed  for  the 
top  and  bottom  slabs. 

e,.  Notes  explaining  the  sign  convention  are  printed  at  the  end  of 
the  output. 

No  summary  sheet  is  made  for  analysis  runs.  In  the  design  phase 
after  all  the  cases  have  been  processed  to  determine  the  most  economic 
thicknesses,  all  cases  are  recomputed  to  obtain  the  moments,  thrusts, 
shears,  etc.,  for  case  output.  During  this  recomputation  the  summary 
sheet  information  is  extracted.  The  summary  sheet  contains: 

a, .  The  first  line,  giving  the  culvert  title  and  date  of 

computation. 

b.  Thicknesses  of  the  top,  side,  and  bottom  numbers. 

c_.  The  maximum  area  of  steel  in  square  inches  required  to  rein- 
force the  inside  and  outside  faces  at  the  seven  computation 
points. 

d.  Range  of  contraflexure  in  the  top  and  bottom  slabs. 

e_.  Cross-sectional  area  of  concrete. 

£.  Diagram  to  depict  the  computation  points. 

Program  number  for  reference. 

This  output  description  is  for  normal  output.  Should  the  user 
have  an  interest  in  the  computation  cycles  that  are  made  or  should  he 
want  to  do  some  debugging,  he  must  put  a C in  the  constant  card,  as 
explained  in  the  input  section.  This  type  of  entry  will  increase  the 
amount  of  output  about  four  times.  The  greater  part  of  the  debugging 
output  is  self-explanatory.  If  the  user  elects  to  run  an  analysis  for 
specific  thicknesses,  no  design  will  be  performed,  nor  will  a summary 
sheet  be  printed. 


31 


Box  culverts 

The  Vicksburg  program  appears  to  be  well  written  and  meets  all  the 
requirements  of  a good  program.  With  a minimum  of  input,  the  program 
produces  the  necessary  information  to  design  a single-cell  box  culvert 
without  haunches. 

Conduits 

None  of  the  conduit  programs  completely  fulfills  the  requirements 
of  a good  design  program.  They  will  all  analyze  a section,  and  the 
Omaha  and  Louisville  versions  compute  the  amount  of  steel  required;  but 
the  concrete  cross  sections  remain  as  they  were  input.  The  Omaha  and 
Louisville  programs  use  the  Working  Stress  Design  method,  while  the 
Southwestern  Division  program  uses  the  Ultimate  Strength  Design  method. 

In  each  case  this  means  the  results  must  be  checked  manually  by  the 
alternate  method.  The  following  is  a summary  of  the  programs. 

Omaha.  The  Omaha  version  is  probably  the  best  of  the  existing 
programs.  It  is  the  only  one  that  allows  multiple  loading  cases  in  a 
single  run.  The  design  mode  is  limited  because  the  concrete  wall 
thicknesses  are  not  sized  by  the  computer.  The  program  does  not  have 
an  analysis  mode. 

Southwestern  Division.  This  is  the  only  program  with  a true  de- 
sign mode,  but  is  is  almost  useless  because  the  wall  thicknesses  are 
sized  by  the  ACI  code  instead  of  the  EM  formulas.  If  the  program  was 
revised  to  size  the  walls  by  the  diagonal  tension  formula  given  in  the 
EM,  it  would  be  a very  good  program.  In  its  present  form  only  the 
analysis  mode  should  be  used.  The  geometry  routines  in  this  program 
are  very  good. 

Louisville.  This  program  was  written  for  a variable  section,  and 
the  designer  is  required  to  input  all  the  voussoir  coordinates  for  the 
concrete  cross  section.  If  a circular  shape  of  uniform  thickness  is  to 
be  analyzed,  one  of  the  other  programs  would  be  easier  to  use.  It  seems 
that  the  curved  beam  correction  factor  is  applied  only  to  the  concrete 
stress  and  not  to  the  steel.  The  program  does  not  have  an  analysis  mode. 
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New  Conduits  program.  If  the  best  parts  of  the  existing  programs 
were  combined.,  a very  useful  program  could  be  developed.  The  Geometry, 
the  Load  Determination,  and  Moment,  Thrust,  and  Shear  phases  of  the 
Southwestern  Division  program  combined  with  the  Design  phase  of  the 
Omaha  or  Louisville  programs  would  be  the  best  combinations.  The  Ulti- 
mate Strength  phase  could  be  used  as  an  investigation  routine  after  the 
section  was  designed  by  the  Working  Stress  method.  When  designing  one 
of  the  conduits  of  constant  thickness,  the  computer  should  increase  the 
wall  thickness  if  either  the  allowable  concrete  design  stresses  or  the 
maximum  area  of  steel  input  is  exceeded.  An  analysis  mode  should  be 
provided  to  investigate  a section  with  predetermined  wall  thickness  and 
area  of  steel.  The  following  is  a suggestion  of  how  the  input  and  out- 


put  should  appear. 

Card 

Format 

Column 

Description 

Specification 

CARD  NO.  1 - Header  information.  The  first  TO  columns  of  this 
card  are  used  for  problem  identification.  Information  contained 
in  the  card  will  be  read  and  listed  at  the  beginning  of  the  print- 
out. The  last  10  columns,  if  used,  will  be  ignored  by  the  com- 
puter. In  the  IBM  1620  version  of  the  program,  this  record  is 
entered  under  a Hollerith  specification.  In  the  GE  225  version, 
it  is  entered  under  a series  of  "A"  specifications  (11A6,  A^+). 

CARD  No.  2 - Control  card. 


10 

1 if  analysis,  zero  if  design 

110 

20 

Number  of  cases  to  be  run 

110 

30 

Print  control 

110 

CARD  No.  3 - Geometry  data.  The  second  record  of  input  data  con- 
tains the  following  eight  values. 

1-5  Number  of  voussoirs  (NV)*  15 

6-10  Number  of  divisions  of  tangent  length  in  15 

oblong  section  (NDIV)** 


* For  circular  conduit,  NV  = 20  ; for  horseshoe  conduit,  NV  = 21  ; 
for  oblong  conduit,  NV  = 20  + NDIV  ; for  nonstandard  section, 

NV  = any  number  not  exceeding  30. 

**  For  any  type  of  section  other  than  oblong,  NDIV  = 0 ; (or  may  be 
left  blank  if  data  is  to  be  processed  by  GE  225). 
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Card 

Column 


Format 

Specification 


Description 


11-15 

Control  parameter  denoting  type  of 
standard  section  (M)* 

15 

16-20 

Control  parameter  denoting  a variable 
section  (MM)** 

15 

21-30 

Thickness  of  section  at  the  crown,  in  feet 

(TC) 

F10.0 

31-1+0 

Inside  radius  of  standard  section,  in  feet 

(RI)  + 

F10.0 

1+1-50 

Tangent  length  of  oblong  section,  in 
feet  ( TANL ) tt 

F10.0 

51-60 

Invert  elevation,  in  feet  (ELINV), 
positive  only 

F10.0 

CARD  Mo.  k - Soils  data.  The  third  record  consists  of  the  fol- 
lowing soils  data  elements. 

1-10  Coefficient  of  vertical  earth  pressure  (VK) 

11-20  Coefficient  of  horizontal  earth  pressure  (HK) 

21-30  Soil  elevation,  in  feet  (SELEV)+ 

31-1+0  Water  elevation,  in  feet  (WELEV)+ 

1+1-50  Rock  line  elevation,  in  feet  (RELEV)+ 

51-60  Unit  weight  of  moist  soil,  in  kef  (WEM) 

61-70  Buoyant  unit  weight  of  saturated  soil,  in 
kef  (WES) 

71-80  Triangular  horizontal  loading,  in  ksf  (THINT) 

CARD  MO.  5 - Design  criteria. 

1-10  Interior  radius  used  in  computing  curved  heam  factor 

11-20  Minimum  cover  of  reinforcement  (blank  causes 
program  to  insert  h inches) 


* For  circular  conduit,  M = 1 ; for  horseshoe  conduit,  M = +1  ; 
for  any  other  type  of  section,  M = 0 (may  be  left  blank  if  data  is 
to  be  processed  by  GE  225). 

**  For  variable  section,  MM  = any  number  other  than  zero;  for  standard 
section,  MM  = 0 (or  blank  if  data  to  be  processed  by  GE  225). 

t For  a nonstandard  section,  RI  = 0 (may  be  left  blank  if  data  is  to 
processed  by  GE  225). 

ft  For  anything  other  than  an  oblong  section,  TANL  = 0 (may  be  left 
blank  if  data  is  to  be  processed  by  GE  225). 

+ All  elevations  must  be  entered  as  positive  values. 


F10.0 

F10.0 

F10.0 

F10.0 

F10.0 

F10.0 

F10.0 

F10.0 
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Card 

Column 

21-30 

31-1*0 

Ul-50 

51-60 

61-70 

71-80 


Format 

Description  Specification 

Minimun  area  of  steel  (blank  causes  program 
to  insert  0.6  sq  in.) 

Maximum  area  of  steel  (blank  causes  program 
to  insert  3.12  sq  in.) 

fy  (blank  causes  program  to  insert  1*0,000  psi) 

f(  (blank  causes  program  to  insert  1*,000  psi) 

Index  of  compressive  strength  of  concrete 
(blank  causes  program  to  insert  0.1*5) 

Required  factor  of  safety  in  Ultimate  Strength 
(blank  causes  program  to  insert  1.8) 


V0USS0IR  DEFINITION  DATA.  The  following  cards  of  the  input  deck 
are  used  to  define  the  shape  of  a variable-section  conduit  and 
are  therefore  not  required  for  the  standard  conduits  of  constant 
thickness.  There  should  be  one  card  for  each  voussoir  of  the 
variable  section,  and  each  card  contains  four  data  elements  as 
follows . 


1-10 

X-coordinate  of  voussoir  extrados, 
in  feet  (X_^) 

F10.0 

11-20 

Y-coordinate  of  voussoir  extrados, 
in  feet  (Y^) 

F10.0 

21-30 

X-coordinates  of  voussoir  intrados, 
in  feet  (Z^) 

F10.0 

31-1*0 

Y-coordinate  of  voussoir  intrados, 
in  feet  (W^) 

F10.0 

Note: 

All  voussoir  coordinates  have  their  origin 
must  be  entered  as  positive  values. 

REMAINING  CARDS.  Card  No.  1*  is  repeated  for 

at  crown  extrados 
each  loading  casi 

Output 

Voussoir  coordinates.  The  coordinates  are  always  printed  on  the 
first  page  of  the  printout. 

Basic  output.  If  a zero  or  a blank  is  input  in  card  column  30  of 
the  print  control  card  (Card  No.  2),  the  following  items  will  be  printed 
for  each  loading  case  (one  case  per  page). 

Input  items  in  the  same  sequence  in  which  they  are  entered 
Voussoir  number  (beginning  at  zero) 

Voussoir  thickness  (inches) 

Shear  (kips) 
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Thrust  (kips) 

Moment  ( kip-f t ) 

Eccentricity  (inches) 

This  would  give  a fairly  clean  printout  and  would  allow  the  de- 
signer to  decide  where  the  changes  in  reinforcing  steel  should  be  made. 

Extended  output.  If  additional  output  is  desired,  punch  a ”1”  in 
column  30  if  the  print  control  card  and  the  following  items  will  be 
printed  (one  case  per  page). 

All  the  items  listed  in  the  Basic  Output  plus  the  following: 
Diagonal  tension  (ksi) 

Area  of  steel  (sq  in.) 

Concrete  stress  (ksi) 

Steel  stress  (ksi) 

Factor  of  safety  in  flexure 
Curved  beam  correction  factor 
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by- 

Carney  M.  Terzian* 

Scope 

This  report  is  an  evaluation  of  the  computer-aided  design  of  con- 
crete iaulticell  conduits.  These  conduits,  typically,  are  used  in  river 
diversion  through  urban  areas  and  in  earthen  dam  embankments  under 
high  fills.  I will  discuss  both  single  purpose  and  general  purpose 
computer  programs. 


Purpose 

The  intent  of  this  evaluation  is  to  aid  the  structural  engineer 
in  facilitating  his  work  while  achieving  the  optimum  analysis  and  de- 
sign. I hope  the  programs  available  and  discussed  here  will  advance 
the  program  of  conduit  analysis  and  design,  along  with  conserving  engi- 
neer time  and  effort. 


Design  Philosophy 


I 


Methodology  is  the  prime  concern-,  two  divisions  are  apparent  in 
the  evaluation  and  will  be  discussed  separately  and  summarized  in  the 
conclusion.  A common  aspect  between  single  and  general  p\irpose  pro- 
grams is  the  criteria  for  design  contained  in  EM  1110-2-2902  Conduits , 
Culverts  and  Pipes  (3  March  1969).  If  we  assume  that  the  EM  repre- 
sents the  state  of  the  art,  then,  presumably,  many  conflicts  are  re- 
solved. The  dissimilarity  between  single  purpose  and  general  purpose 
programs  has  been  considerably  reduced  and  one  is  as  good  as  another. 
But,  consider  the  state  of  the  art.  Three  important  references  in  the 
field  are: 


* Chief,  Structural  Section,  hew  England  Division. 


37 


a.  "Finite  Element  Analysis  of  Port  Allen  and  Old  River  Locks”  by 
Clough  & Duncan,  Sep  1969,  University  of  California,  Berkeley. 

b_.  "Thick  Walled  Multiple  Opening  Reinforced  Concrete  Conduits" 
by  Ryan,  Salem,  Gamble  & Mohraz,  Dec  1972,  University  of 
Illinois. 

c_.  "Non-Linear  Analysis  of  Planar  Reinforced  Concrete  Structures" 
by  Salem  & Mohraz,  Jul  197^,  University  of  Illinois. 

The  findings  of  these  Corps  sponsored  researches  have  not  yet  been 
incorporated  by  the  subject  EM.  For  example.  Plates  3 and  4,  note  3,  of 
the  EM  specifies  the  assumption  of  uniform  foundation  pressures;  Ref- 
erence (a)  suggest  yielding  foundation  and  its  Appendix  C contains  a 
finite  element  computer  program  for  Soil-Structure  Interaction.  Ref- 
erences (b)  and  (c)  include  information  developed  on  the  use  of  the 
finite-element  method  of  analysis,  taking  into  account  the  progression 
of  cracking  with  increasing  load  and  the  nonlinear  stress-strain  re- 
sponse of  the  materials.  Obviously  the  researchers  are  using  the  finite- 
element  method,  the  basis  for  general  purpose  programs.  The  subject  EM 
should  be  revised  to  recognize  these  advances  and  others  as  they  are 
developed. 


Evaluations  of  Single  and  General  Purpose  Programs 
Single  purpose  programs 

The  rating  procedure  for  single  purpose  programs  is  shown  in 
Table  1 . 

Rating  of  single  purpose  programs 

The  single  purpose  programs  were  evaluated  using  the  rating  pro- 
cedure shown  in  Table  1.  Checking  the  fifteen  elements  considered 


necessary  for  optimum 
normalized  grading: 

analysis  and  design  resulted 

in  the  following 

Grade 

Program 

District 

A+ 

713-G1-M1070 

Albuquerque 

A 

13-G1-AU07 

St.  Louis 

B+ 

713-G2-L2-005 

San  Francisco 

D+ 

713-x6-L2-22A 

Sacramento 

Discussion  of  ratings 


In  general  the  approach  to  analysis  by  moment  distribution  was 
similar  for  all  programs  rated  (except  L2-005);  the  difference  occurs 
in  the  design  function.  Program  No.  713-x6-L2-22A,  Sacramento,  as 
transmitted  provided  no  design  of  elements,  resulting  in  a low  rating. 
Genral  purpose  programs 

The  rating  procedure  for  general  purpose  programs  is  shown  in 
Table  2. 

Ratings  of  general  purpose  programs 

The  general  purpose  programs  were  evaluated  using  the  rating  pro- 
cedure given  in  Table  2.  Checking  the  fifteen  elements  considered 
necessary  for  optimum  analysis  and  design  resulted  in  the  following 


normalized  gradings: 

Grade 

Program 

District 

A+ 

STRUDL/MCAUTO 

Proprietary 

B 

SAP  IV 

University  of  California 
Berkeley 

C 

EFFRAM (713-D0F7110) 

New  England  Division 

C- 

Wilson  2D(713-G2-L3-002) 

San  Francisco 

D 

G-FRAME ( 713-F3-A1- 030 ) 

Memphi s 

Discussion 

of  ratings 

All  general  purpose  programs  used  the  forced  displacement  finite- 
element  method.  The  primary  differences  in  program  capabilities  is 
that  MCAUTO  provides  3-dimensional  and  dynamic  analysis,  interactive 
graphic  (CRT),  and  design  while  SAP  IV  did  not  contain  the  design 
feature.  The  remaining  three  programs  were  rated  low  since  only  two- 
dimensional  analysis  is  performed,  with  no  dynamic  analysis  and  no  in- 
teractive graphics  or  design  capabilities  as  yet.  It  is  noted  that 
STRUDL/MCAUTO  with  FASTDRAW  is  a total  system  while  the  others  are 
programs  only.  The  low  ratings  are  not  completely  fair  since  EFFRAM, 
Wilson  2D,  and  G-FRAME  provide  excellent  analytical  and  simulation  pro- 
cedures not  found  in  the  single  purpose  programs.  They  could  be  rated 
higher  if  compared  to  the  single  purpose  programs. 
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Two  views  of  the  computer  programs  will  be  made,  first  by  group- 
ings, then  followed  by  the  best  points  of  each  program. 


Four  general  areas  of  concern  were  defined  in  the  comparison  of 
the  nine  programs;  these  are  simulation,  execution,  design,  and  results. 
Specifically,  items  shown  in  the  rating  procedure  are  grouped  as  follows: 


a.. 

Simulation  - Items  1,  5,  6, 

7,  10,  12,  lit,  & 15 

b. 

Execution  - Items  3,  1+,  9, 

& 11+ 

c_. 

Design  - Item  8 

d. 

Results  - Items  2,  11,  13, 

& ll+ 

The  top  3 of  9 programs  are  rated,  based  upon  these  four  groupings, 
as  follows: 

a.  Simulation  - MCAUTO,  SAP  IV,  & EFFRAM 

l).  Execution  - MCAUTO,  SAP  IV,  & EFFRAM 

c.  Design  - MCAUTO,  M1070,  & Al+07 

d.  Results  - MCAUTO,  SAP  IV,  & M1070 

It  has  been  found  that  each  of  the  programs  studied  has  its  own 
unique  qualities  that  set  it  apart  in  one  category  or  another  from  all 
the  others.  We  feel  that  it  would  be  useful  to  all  of  you,  in  your 
future  program  development,  to  keep  these  following  ideas  in  mind. 

Foundation  springs,  capable  of  taking  tension  or  compression, 
should  be  provided  with  the  option  of  eliminating  any  tension  that 
occurs  in  order  to  accurately  simulate  any  unique  environmental  con- 
ditions that  may  arise.  This  ability  is  present  in  EFFRAM  by  NED. 

Along  the  same  reasoning,  torsional  or  moment  springs  should  be  pos- 
sible, as  in  SAP  IV  or  Wilson  2D.  A yielding  foundation,  we  feel,  is 
an  essential  element  in  any  totally  comprehensive  structural  analysis 
program . 

The  program  should  be  capable  of  accepting  as  many  types  of  load- 
ing as  possible  without  the  need  to  compute  FEM's  or  resultant  applied 
moments  in  place  of  the  actual  loading  itself.  EFFRAM  gives  the  engi- 
neer the  ability  to  apply  surcharge,  strip,  and  point  loads  at  varying 
distance  from  a vertical  wall  as  well  as  uniform,  concentrated,  and 

1+0 
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distributed  loadings.  FEM's  are  provided  by  the  program  rather  than  the 
engineer,  therefore  reducing  the  chance  of  error  and  the  time  needed  to 
set  up  and  check  input  data. 

Program  input  and  output  should  be  arranged  so  that  an  engineer 
without  a great  working  knowledge  of  the  program  itself  can  set  up  in- 
put data  or,  on  the  other  hand,  be  able  to  readily  interpret  output.  It 
is  therfore  advised  that  abreviations,  codes,  or  symbols  be  avoided. 
Labels  should  be  given  so  as  to  enable  all  engineers  to  understand  the 
type  of  information  being  provided.  Results  can  therefore  be  easily 
circulated  throughout  the  country  without  the  need  for  time  consuming 
explanations  which  are  usually  not  fully  understood.  Standard  engineer- 
ing notation  should  be  used.  G-FRAME  is  a good  example  of  this.  It  is 
simply  arranged  and  set  up  for  everyone. 

Programs  that  have  design  capabilities  should  contain  as  many 
"default"  characteristics  as  possible  so  that  when  standard  engineering 
parameters  are  used  there  will  be  no  need  to  input  them.  STRUDL  provides 
this  excellent  capability;  it  can  also  be  seen  to  a limited  degree  in 
the  programs  "Concrete  Box  Culvert  Frame  Analysis  and  Design"  (T13-G1- 
M1070)  and  "Multi-Cell  Box  Culvert"  (13-G1-AU07 ) . 

SAP  IV  is  unique;  it  will  handle  dynamic  as  well  as  static  analy- 
sis. This  feature  is  rarely  seen  in  a single  program  and  expands  its 
potential  use  considerably. 

Another  very  useful  program  aid  is  member,  node,  and  load  genera- 
tion. This  capability  is  of  particular  value  in  programs  using  finite- 
element  techniques  where  a great  many  node  coordinates  and  member  prop- 
erties and  incidences  must  be  input  to  simulate  the  design  problem. 
Through  this  generation  technique,  the  engineer  is  able  to  save  con- 
siderable amounts  of  time  setting  up  and  inputting  his  design  data. 

SAP  IV,  Wilson  2D,  and  EFFRAM  exhibit  this  very  useful  characteristic. 

The  single-  or  multiple-culvert  analysis  by  Los  Angeles  District 
has  a unique  idea  for  outputting  program  data.  Conduit  moments  are 
printed  in  a diagrammatic  shape  of  the  structure  under  analysis.  This 
method  of  formatting  output  provides  the  engineer  with  his  required 
moments,  showing  their  respective  locations,  while  at  the  same  time 
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depicting  the  shape  of  the  structure  so  that  the  engineer  can  readily 
keep  it  in  mind.  All  this  is  done  in  one  simple  step. 

With  regard  to  any  design  procedures  in  addition  to  a pure  analy- 
sis, STRUDL  is  a prime  example  of  the  capabilities  that  could  be  pro- 
vided. The  option  of  designing  in  either  concrete  or  steel  by  WSD  or 
USD  is  present  along  with  numerous  code  checking  procedures.  Member 
property  tables  are  also  provided.  The  engineer  need  only  call  out  the 
structural  shape  he  desires  without  the  need  of  determining  member 
properties,  such  as  area,  section  modules,  or  moment  of  inertia,  to 
name  just  a few. 

Graphic  plotting  and  interactive  graphics  are  other  relatively 
new  areas  which  many  engineers  have  yet  to  explore.  With  the  aid  of 
those  tools,  the  engineer  acquires  the  ability  to  put  his  resultant 
design  into  an  advanced  analytical  perspective.  He  has  the  ability  to 
display  his  original  model,  graphically  providing  him  with  an  efficient 
and  economical  way  of  checking  his  input  data.  He  can  ask  for  plots  of 
moments,  shears,  or  influence  lines  about  any  axis  or  display  a two- 
or  three-dimensional  deflected  shape  of  his  structure  under  a particular 
loading  condition.  Through  interactive  graphics  he  can  model  a struc- 
ture with  the  aid  of  a tablet  and  input  pen.  He  can  change  the  shape 
of  a structure  instantly  and  ask  for  a reanalysis  without  having  to 
rearrange  members  and  coordinates  manually.  The  possibilities  are 
limited  only  by  the  imagination  of  the  engineer  himself.  McDonnel 
Douglas  Automation  Company  provides  this  service  with  its  propietary 
STRUDL  package.  The  program,  "Concrete  Box  Culvert  Frame  Analysis  and 
Design,"  attempts  to  achieve  some  of  the  advantages  graphics  may  provide 
by  modeling  the  structural  shape  under  consideration  in  its  output 
format . 

Finally,  one  last  characteristic  which  all  the  programs  exhibit 
should  be  commented  on.  FORTRAN  has  been  chosen  as  the  programming 
language  in  all  cases.  It  really  is  the  only  appropriate  language  to 
use  when  dealing  with  engineering  or  scientific  principles.  It  has 
been  found,  though,  that  some  programs  in  the  WES  library  are  present 
in  whole  or  in  part  in  a language  peculiar  to  the  machine  that  they 
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were  written  on.  In  order  to  achieve  the  ends  of  a Corps-wide  library, 
programs  must  be  written  so  that  they  are  freely  accessible  to  all 
potential  users  without  the  need  to  perform  large  conversion  procedures 
that  will  put  the  program  in  a "ready  to  run"  stage.  Those  in  charge 
of  the  Corps  computer  systems  should  be  able  to  provide  the  engineer 
who  is  constantly  engaged  in  compter-aided  design  with  the  procedures 
necessary  to  run  a program  on  a particular  computing  system.  At  present, 
this  information  has  not  been  found  to  be  readily  available. 

Notes  on  Specialty  Session  C 

The  session  clarified  the  importance  of  simulation  by  stressing 
the  use  of  the  finite-element  method,  soil-interaction,  and  computer 
graphics.  The  Design  phases  of  most  programs  are  limited  to  the  extent 
that  the  design  may  be  for  USD  and  not  WSD  or  vice  versa.  None  of  the 
programs  have  Automatic  Design  based  on  Corps  criteria  (EM  1110-1-2101). 
Improvement  in  the  Design  phase  is  imperative. 

The  session  also  indicated  the  predominant  use  of  time-shared 
terminals  with  acoustic  pickups  to  WES  or  leased  computer  services 
such  as  CSC  and  MCAUTO. 


TUMNELS 


by 

Robert  J.  Smith* 
Introduction 


This  paper  presents  the  "State-of-the-Corps-Art"  for  computer- 
aided  design  of  tunnels.  A tunnel  may  be  defined  as,  "an  underground 
passage  made  without  removing  the  overlying  rock  or  soil."  Most  Corps 
of  Engineers  tunnels  are  either  water  conveyance  or  transportation 
tunnels  driven  through  rock.  Therefore,  this  paper  is  limited  to  ways 
in  which  the  computer  can  be  used  to  aid  the  structural  engineer  in  the 
design  of  such  tunnels. 

The  structural  engineer  must  work  closely  with  the  geologist  and / 
or  soils  engineer  and  hydraulic  engineer.  The  structural  engineer's 
responsibility  is  to  design  the  temporary  and  permanent  supports  neces- 
sary to  stabilize  the  tunnel.  Tunnel  linings  are  designed  using  a com- 
bination of  theory,  intuition,  and  experience.  As  a result,  tunnel  lin- 
ing design  is  more  of  an  art  than  a science.  Thus  coordination  with 
other  disciplines  is  very  important  when  designing  tunnel  support 
systems. 


Types  of  Support 

As  stated  above,  the  structural  engineer  is  responsible  for  the 
design  of  temporary  and  permanent  supports.  The  temporary  support  may 
consist  of  any  one  or  combination  of  the  following:  (a)  steel  sets, 

(b)  rock  bolts,  (c)  spiling,  or  (d)  shotcrete.  The  permanent  support 
could  be  any  one  or  combination  of  the  above,  but  for  a water  convey- 
ance tunnel  it  would  probably  be  concrete.  If  the  tunnel  lining  is  to 
withstand  internal  water  pressure,  a steel  lining  inside  the  concrete 
lining  will  probably  be  required. 


* Structural  Engineer,  OCE. 


Programs  Available 


Now  that  we  have  identified  the  structural  engineer's  responsibil- 
ities and  types  of  supports  used  to  stabilize  a tunnel,  let  us  look  at 
the  computer  programs  that  are  available. 

Computer  program  for  steel  sets 

Steel  sets  may  be  designed  using  the  program  presented  in  Water- 
ways Experiment  Station  (WES)  Technical  Report  C-73-2,  "Computer  Study  of 

Steel  Tunnel  Supports."'*’  This  program  is  written  using  the  philosophy 

2 

presented  by  R.  V.  Proctor  and  T.  L.  White.  This  method  of  design  is 

R 

also  presented  m the  draft  of  EM  1110-2-2901.  The  Proctor  and  White 

method  is  basically  a graphical  solution.  The  steel  set  is  assumed  to 

be  pinned  at  each  blocking  point  with  each  blocking  point  carrying  its 

proportionate  share  of  the  rock  load.  The  rock  load  is  determined  by 

using  the  Terzaghi  rock  load  concepts.  The  next  step  is  to  construct 

a force  polygon.  From  this,  stresses  in  the  ribs  can  be  determined. 

The  method  is  illustrated  in  the  two  references  mentioned  above. 

The  computer  program  "Tunnel"  uses  information  about  geometry, 

loads,  location  of  blocking  points,  and  member  properties  to  obtain 

shears,  moments,  thrusts,  and  displacements.  The  program  uses  the 

stiffness  matrix  method  for  analyzing  the  steel  sets.  The  program  is 

designed  to  handle  all  tunnel  shapes  and  is  written  in  a manner  which 

will  allow  the  user  to  have  no  knowledge  of  computer  programming  and 

very  little  knowledge  of  his  operating  systems. 

There  are  several  disadvantages  to  the  graphical  solution  pre- 

2 

sented  by  Proctor  and  White.  The  analysis  yields  no  information  on 
deflections  within  the  structure.  The  method  does  not  account  for 
yielding  (elastic)  behavior  of  the  steel  sets.  Both  of  the  above  con- 
ditions affect  the  loads  acting  upon  the  structure.  Furthermore,  the 
graphical  method  is  tedious  and  time-consuming  and  causes  the  structural 
engineer  to  avoid  consideration  of  several  alternative  designs. 

The  computer  program  is  designed  to  alleviate  many  of  the  dis- 
advantages of  the  above  method.  Although  the  loads  are  treated  as 
being  independent  of  the  distortions  in  the  system,  the  program  is 
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flexible  enough  to  accept  any  type  of  loading.  Not  only  is  information 
obtained  on  the  deflections  of  the  structure,  but  member  moments  and 
shears  as  well  as  member  thrusts  are  obtained  directly.  The  speed  with 
which  the  computations  are  performed  allows  the  designer  freedom  to 
make  and  choose  from  several  designs. 

"Tunnel"  is  written  in  FORTRAN  IV  for  a time-sharing  computer 
system.  The  report  contains  an  auxiliary  program  called  "Horshu"  which 
is  used  to  automatically  generate  the  required  information  for  "Tunnel." 
The  WES  report  consists  of:  an  introductory  section  that  presents  cur- 

rent steel  set  tunnel  support  design  procedures;  a description  of  the 
main  program,  "Tunnel";  and  a description  of  the  auxiliary  program, 
"Horshu."  The  mathematical  basis  for  "Tunnel"  is  discussed  in  Appendix 
A.  Appendix  B contains  a blank  sample  data  input  form.  Appendix  C is 
a listing  of  the  program  "Tunnel"  and  Appendix  D is  a listing  of  the 
program  "Horshu." 

The  programs  have  not  been  used  for  design;  however,  the  Corps 
has  39  tunnels  in  the  planning  stages,  so  perhaps  they  can  be  used  in 
the  future. 

Computer  programs  for  concrete  linings 

Responses  to  the  questionnaire  revealed  that  four  different  com- 
puter programs  have  been  used  for  design  of  concrete  tunnel  linings. 

The  computer  program  numbers  and  district  offices  that  have  used  them 
are  as  follows: 

District  Office  Program  Number 


Albuquerque 

Huntington 

Portland 

Sacramento 


713-G1-M0090 

713-K2-H1325 

513-06-00321 

713-61-L202A 


The  philosophy  used  in  these  programs  was  the  same  as  discussed 
earlier  in  the  conduit  part  of  this  specialty  session  (or  in  the  papers 
on  conduits  written  for  this  specialty  session).  Therefore,  it  will 
not  be  discussed  here;  but  some  of  the  loading  assumptions  that  must  be 
considered  when  using  one  of  these  programs  for  designing  concrete 
tunnel  linings  are  as  follows:  (a)  resist  hydrostatic  water  load 

acting  alone,  (b)  resist  rock  load  acting  alone,  (c)  resist  a 
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combination  of  the  water  and  rock  loads,  and  (d)  resist  the  grouting 
pressures. 

The  above-listed  loading  conditions  are  not  to  be  construed  as  all 
inclusive  but  rather  are  given  so  the  reader  can  have  an  idea  of  the 
minimum  number  of  runs  to  make  when  designing  concrete  tunnel  linings. 

Computer  program 

for  steel  tunnel  liners 

Sacramento  District  has  used  program  number  713-G1-L202C,  "Liner," 
to  design  the  cylindrical  steel  tunnel  liners  used  for  the  New  Melones 
f Project.  For  external  pressure  the  shell  is  assumed  to  be  confined  in 

a rigid  cavity;  stiffener  rings  are  incorporated  as  needed.  For  inter- 
nal pressure,  division  of  load  between  liner  and  rock  is  based  upon 

OCE  criteria  established  for  this  project.  Plane  strain  is  assumed. 

k 

Draft  of  EM  1110-2-2907  contains  a discussion  on  the  considerations  to 
be  given  for  this  type  of  design.  The  program  can  be  altered  to  sat- 
isfy the  design  requirements  of  a specific  tunnel. 

Finite-element  method 

The  philosophy  of  this  type  of  program  will  not  be  discussed  here 
since  it  is  documented  in  many  other  papers  and  books.  What  will  be 
discussed  are  considerations  to  be  made  when  using  this  type  of  analy- 
sis. An  approach  that  can  be  taken  is  as  shown  in  the  following  outline: 

I.  Isolate  an  appropriate  geostructural  block  containing  the 
opening  or  excavation  from  the  rock  mass  for  structural 
analysis. 

II.  Determine  the  geology  of  the  geostructural  block. 

III.  Determine  forces  (or  stresses)  acting  upon  or  within  the 
geostructural  block. 

A.  Surface  forces 

1.  External  gravity  (including  lateral  rock  pressure) 

2.  Tectonic 

3.  Groundwater 
U . Earthquake 
5 . Man-made 

B.  Body  forces  within  the  model 
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1.  Internal  gravity 

2.  Residual 

3.  Groundwater 

4 . Earthquake 

5.  Man-made 

C.  Due  to  opening  or  cut 

D.  Due  to  varying  stiffness  of  the  rock 

E.  Due  to  relaxation  or  creep 

IV.  Determine  the  mechanical  rock  properties. 

V.  Develop  the  grid. 

VI.  Determine  if  failure  conditions  (fracture,  slippage,  excessive 

deflection)  exist. 

A.  No  failure  conditions  - analysis  is  complete.  No  struc- 
tural support  required. 

B.  Failure  conditions  - some  type  of  structural  support 
required. 

VII.  Determine  structural  support  required  to  stabilize  opening. 

The  above  outline  was  extracted  from  a paper  written  by 
J.  S.  Dodd.5 

Another  discussion  on  the  finite-element  method  of  analysis  can 
be  found  in  a report  prepared  by  J.  J.  K.  Daemen.^  This  report  dis- 
cusses development  of  ground  reaction  curves,  stiff  and  soft  support- 
ing methods,  and  limiting  convergence  of  the  tunnel  walls.  A computer 
program  called  Tunsup  is  presented  in  the  report.  This  program  goes 
beyond  the  analysis  presented  above  by  considering  things  such  as 
ground- support  interaction,  incremental  loading  or  unloading  (equivalent 
mining),  and  strain-softening  of  the  rock  around  the  tunnel  periphery. 
The  program  is  written  in  FORTRAN  IV.  For  the  version  presented, 

135>000  central  memory  words  are  required  on  a Control  Data  Corporation 
CYBER  74  computer  with  KRONOS  operating  system  and  RUN  compiler.  Run- 
ning times  vary  greatly  with  problem  type.  Computer  times  for  an  ex- 
ample given  in  the  report  ranged  from  2 up  to  l8  minutes. 

The  program  has  not  been  used  for  design  but  may  have  a place  in 
the  design  of  the  more  complex  tunnels.  Use  of  "Tunsup"  is  presented 
in  Appendix  C of  the  MRD  report. 
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Programs  Needed 

Rock  bolts,  spiling,  and  shotcrete  have  not  been  discussed  thus 
far.  These  support  systems  could  be  modeled  using  the  finite-element 
method  of  analysis.  However,  the  short  literature  search  that  was  made 
in  preparing  this  paper  did  not  reveal  any  designs  that  had  used  the 
computer  as  a design  aid  for  these  types  of  supports. 

Conclusions 

As  was  stated  in  the  Introduction,  tunnel  lining  design  is  more 
of  an  art  than  a science;  thus,  the  computer  cannot  be  used  to  automate 
the  design  of  tunnel  supports.  The  computer  can,  however,  be  used  to 
assist  in  making  sound  engineering  judgments.  The  size  of  steel  sets, 
thickness  of  concrete  linings,  spacing  and  length  of  rock  bolts,  and 
thickness  of  shotcrete  should  not  vary  too  much  from  the  commonly  used 
"rules  of  thumb."  It  should  be  remembered  that  the  input  is  not  ^xact, 
therefore,  the  temporary  and  permanent  supports  cannot  be  designed 
down  to  the  last  pound  of  steel  or  cubic  yard  of  concrete.  European 
trends  (so  called  New  Austrian  Tunneling  Method)  are  toward  emphasis  on 
observation  and  control  of  lining  displacements  as  the  primary  variable 
in  design.  They  put  less  emphasis  on  design  before  tunneling  starts. 
United  States  and  British  procedures  have  tended  to  be  concerned  with 
loads,  and  advance  designs  to  resist  loads.  The  structural  engineer 
should  remember  that  the  lining  design  can  be  strengthened  (supple- 
mented) with  additional  bolting,  shotcreting,  and/or  jump  sets  when  ob- 
servations dictate.  Both  design  concepts  should  be  kept  in  mind  when 
designing  tunnel  supports. 

Postconference  Comments 


The  following  comments  and/or  recommendations  were  made  by 
participants : 

a.  St.  Louis  District  is  using  the  WES  computer  program  "Tunnel" 
as  a design  aid  on  Meramec  Park  tunnel. 
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The  Bureau  of  Mines  and  Bureau  of  Reclamation  have  computer 
programs  which  may  he  useful. 

Sacramento  District  program  number  713-61-L202A  should  be  de- 
leted from  the  list  of  programs  which  may  be  used  for  concrete 
Dining  design.  This  was  a recommendation  by  Mr.  Haavisto  of 
the  Sacramento  District. 
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